Zeeman effects on the impurity-induced resonances in rf-wave superconductors 
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It is shown how the resonant states induced by a single spinless impurity in a dj.2_y2-wave super- 
conductor evolve under the effect of an applied Zeeman magnetic field. Moreover, it is demonstrated 
that the spin-orbit coupling to the impurity potential can have important and characteristic effects 
on the resonant states and their response to the Zeeman field, especially when the impurity is close 
to the unitary limit. For zero or very small spin-orbit interaction, the resonant states becomes Zee- 
man splitted by the magnetic field while when the spin-orbit coupling is important, new low-lying 
resonances arise which do not show any Zeeman splitting. 
PACS number(s): 74.20.-z, 7f.70.Ej, 74.62.Dh 
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Recent scanning tunneling microscopy (STM) mea- 
surements in BSCCO high-Tc superconductors have 
recorded clear images _pf quasiparticle resonant siates 
around intrinsic defectsa and individual Aucl and ZnEI im- 
purity atoms. A common result of these different mea- 
surements is that, in average, the resonances are observed 
just below the Fermi energy, indicating a resonant energy 
ujQ much smaller than the superconducting energy gap 
A {\ujo\ ^ A/30). For the case of Zn- impurity doping, 
the spatial dependence of the resonant states exhibits a 
large signal at the impurity site with local maxima on 
the second neighbour Cu sites [i. e. along the node-gap 
directions) followed by somewhat weaker peaks along the 
directions of the gap maxima.l3 

Low-lying quasiparticle resonant states have been pre- 
dicted to occur in dx^^y2-wave superconductors doped 
with spinless impurities close to the unitary limit D Ac- 
cording to such model, the resonance at — 1-5 meV for 
Zn-substituted samples implies a scattering parameter 
c = I/7riVFyii„p of about O.ljj where l^mp is the s-wave 
impurity potential and Np is the normal-state density of 
states for each spin species at the Fermi level. Similar 
values of c are estimated for Au impurities and generic 
intrinsic defectsEfl Furthermore, the observed power-law 
decay G{r) ~ ^/r^'^^ of the angle-averaged differential 
conductance G{r) at large distances r from the impurity 
sitetl is in very good accord with G(r) ^ l/r^, predicted 
in Ref. I 

Despite of the agreements between theory and experi- 
ments, the spatial dependence of the resonant states re- 
ported in Ref. |^ is quite ad odds with that expected by 
strong spinless impurities. These, in fact, would gener- 
ate, in addition to the peak at the impurity site, a four- 
fold symmetric signal with maxima along the directions 
where the gap is fully opened .qIj 

The images recorded by the STM measurements ap- 
pear therefore rotated by 7r/4 with respect to those re- 
sulting from a spinless quasi-unitary impurity potential. 
Recently, it has been proposed that a blocking effect of 
the BiO and SrO layers between the tunneling tip and 
the Cu02 layer would actually give rise to. the same spa- 
tial dependence recorder in experiments. □ On the other 



hand, according to a recent theoretical analysisj^ the ob- 
served spatial dependence would rather arise from a lo- 
cal antiferromagnetic spin re-arrangement induced by a 
nominal zero-spin weak impurity. From this perspective, 
the Zn-atom behaves effectively as a magnetic impurity 
with non-local coupling to the charge carriers leading to 
the X-shaped geometry of the resonant state. Such a 
picture would be consistent also with the presence of un- 
paired S = 1/2 moments in the, vicinity of Zn atoms as 
observed by NMR experiments.^ 

From the above discussion and the contrasting claims 
reported in recent literature, it appears that the problem 
of deciding whether nominal spinless impurities behaves 
as non-magnetic or effectively magnetic scattering cen- 
ters in high- Tc cuprates is still an open issue. However, 
models based on purely non-magnetic impurity poten- 
tials predict resonant states quite seasitive to impurity 
strength and charge carrier doping,Qlj while the picture 
proposed in Ref. ^ has been claimed to yield results much 
more robust. This qualitative difference could therefore 
be used as a tool discriminating between the two pictures 
by, for example, studying the response to some external 
applied perturbation. 

The aim of this paper is two-fold. First, it is shown 
how the resonant states induced by a single spinless im- 
purity in a dx2_y2-wave superconductor evolve under the 
effect of an applied Zeeman magnetic field. Second, it 
is demonstrated that the spin-orbit coupling to the im- 
purity potential becomes especially important when the 
impurity is close to the unitary limit ,113 and can have im- 
portant effects on the resonant states and their response 
to a Zeeman field. According to whether the spin-orbit 
scattering is irrelevant or not, the Zeeman response of the 
resonant state behaves in two distinct ways. For zero or 
very small spin-orbit interaction, the resonances become 
Zeeman splitted and the spatial dependence can change 
from electron-like to hole-like for already quite small val- 
ues of the imposed magnetic field. On the contrary, 
when the spin-orbit coupling to the quasi-unitary impu- 
rity becomes relevant, new low-lying sharp resonances 
arise which do not show Zeeman splitting. In this lat- 
ter case, the spatial dependence at fixed energy can be 
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made to change from electron-like to a novel, spin-orbit 
induced, symmetry for a suitable value of the Zeeman 
field. 

The differential conductance recorded in a STM ex- 
periment is proportional to the local density of states 
(LDOS) N{r, oj) - -(1 A)Im[Gi\(r, r; u;) + G^r, r; c.)], 
where G^'{r,r;uj) is the retarded 4x4 matrix Green's 
function defined in the particle-hole-spin space and r 
is the vector position with respect to the impurity lo- 
cated at the origin. The (11) and (22) components re- 
fer to the two pseudospin states of the quasiparticles. 
For the single impurity case, G^{r,r;uj) is obtained by 
the Fourier transform of GCk^k'^uj) — (5k,k'Go(k, cj) -|- 
Go(k,a;)T(k,k';w)Go(k',cj), where Go{k,uj) is the 
Green's function for the pure ci2^2_j^2-wave supercon- 
ductor and r(k, k';a;) is the T-matrix associated to 
the impurity scattering: r(k, k';a') — F(k, k') -|- 
J2k»V{k,k")Go{k",oj)T{k",'k';uj), where V{k,k') is 
the impurity potential in momentum space. In the fol- 
lowing, an external magnetic field H is assumed to be 
directed parallel to the conducting x-y plane, for ex- 
ample H II X, and the spins are quantized along the 
direction of H. Under the assumption of strong two- 
dimensionality, the coupling of the planar magnetic field 
to the quasiparticle spins becomes predominantj-over 
the coupling to the quasiparticle orbital motion.O In 
the limitips case for which the orbital coupling can be 
neglected Jl3 the Green's function Go(k, uj) is simply given 
by G(7^(k, w) = Lj-€{k)p3-A{k)p2T2 + hp3T3, where e(k) 
is the quasiparticle dispersion, A(k) = A((/>) = Acos(20) 
is the gap function, and h — fi^H is the Zeeman en- 
ergy {hb = Bohr magneton). The Pauli matrices pi 
and Tj = 1,2,3) act on the particle-hole and spin 
subspaces, respectively. For sufficiently low values of 
h/A, the Fulde-Ferrel-Larkin-Ovchinnikov state can be 
ignored, and the effect of H is merely to split the^seudo- 
spin degeneracy of the quasiparticle excitations .Ej~c3 

The impurity atom is assumed to have a simple S- 
function potential: 14nip;(r) — Mmp^(r). According to 
the Elliott- Yafet theory ,E3 the spin-orbit coupling to 14mp 
can be modelled as 14o(r) = {Sg / kp)['Wimp{r) x p] • cr, 
where p = — iV and cr are respectively the momentum 
and spin operators, kp is the Fermi mom en .t um and dg 
is of the order of the shift of the g-factor.tllE^ Here, Sg 
is treated as a free parameter, however not exceeding 
(5^ ~ 0.1 — 0.2, whieh is the expected order of magnitude 
for Cu02 systems.im Since the charge carriers are con- 
fined to move on the x-y plane, only the cr^ component 
of Vsoiic) is nonzero. Hence, in the particle-hole spin sub- 
space, the total impurity potential F(k, k') reduces to 
V(k,k') = + i(55V^mp[k x k']^Ti. 

The k X k' dependence of the spin-orbit contribu- 
tion permits to decouple the T-matrix into two com- 
ponents: r(k, k';(jj) = ri,„p(aj) -I- rso(k, k'; w), where 

71inp('^) = T4mpP3 + Z]k^mp/03Go(k,tj)Timp(w) is the 

usual T-matrix for the scalar potential and: 
Tso(k,k';a;) = iSgVin^p[k x k'J^ri 
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FIG. 1. LDOS without spin-orbit interaction as a function 
of uj/A for c = 0.08 and ft = (solid lines), h/A = 0.02 
(dotted lines), and h/A = 0.04 (dashed lines), (a): LDOS at 
the impurity site r = (0, 0). (b): LDOS along the direction of 
gap maxima r = (4/A:f,0) 
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k X k"].TiGo(k",c^)r,o(k",k';c^), (1) 



is the Tjpatrix for the spin-orbit coupling to the 
impurity.EHl The resulting LDOS is therefore the sum of 
three contributions: N{r,uj) = No{uj) -I- (5iVimp(r, w) -|- 
6Nso{r,uj). From now on, the LDOS contributions are 
given in units of the normal-state DOS summed over the 
two spin directions, 2Np, and particle-hole symmetry is 
assumed. Hence, Nq{uj) — [lmgo{uj+) + lmgo{uj^)]/2, 
where go(^±) = / ^w/[A((/))2 - and w± = 

u} T h, is the LDOS for the pure superconductor, while 
'^iVi„p,(so)(r,a;) = SN+^^^^^^^{r,uj)+SNr^^^^^^^{r,uj) is the 
LDOS contribution induced by the interaction with the 
impurity (spin-orbit) potential. The impurity part of the 
LDOS is just the superposition of the zero-field LDOS 
reported in ReL ^ shifted by ±h: 
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go('',^±)^ ^ /o(r,^± 
5o(t^±) + c go{uj±) ' 
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where 



ffo(r,^^±) 
/o(r,^±) 



2tT [A((/))2 -cj2]1/2 



A(0) 



(2) 
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A crucial effect of the Zeeman magnetic field is that the 
poles arising from the denominators of Eq. (|^) are splitted 
by ft,. In fact, for small values of c and h, the energy 
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FIG. 2. Spatial dependence of the LDOS without 
spin-orbit interaction for c = 0.08 and lu/A — —0.03. (a): 
ft = 0. (b): h/A = 0.02. 



contemporary suppression of the central peak at r = 
and the signals along the gap- node directions. This pat- 
tern is equal to that reported in Fig.|| of Ref. ||, but here 
it has been obtained without reversing the sign of w. As 
can be also inferred from Fig. 1 , higher values of h mod- 
ify the intensity of the signal, but its spatial dependence 
remains equal to that of Figj^b. It is important to stress 
that Fig.|| refers to the spatial dependence on the Cu02 
layer and no blocking effect has been considered. 

Let us consider now under which conditions the re- 
sults of Figs 1 and 2 are modified by the presence 
of spin-orbit coupling to the impurity. The spin-orbit 
T-matrix contribution, Eq.(jl[), is obtained by setting 
Tso(k,k';a;) = i(5.gMmp[k x t(k',w)]^Ti, where t{k,uj) = 
h + iSgVimpEk' k'TiGo(k',cj)[k' xt(k,cj)],y The equa- 
tion for t(k, Lu) is easily solved in terms of its components, 
tx and ty , and after some algebra the resulting spin-orbit 
part of the LDOS reduces to: 



SN^Jr,^) = - -Im{A±H[5e(r,w±)2 +g,(r,c^±) 

+ /c(r,o;±)2+/,(r,o.±)2]} 

- Im{B=^ {cj) [gs (r, uj±)fs (r, uj± ) 

-5c(i-,w±)/c(r,w±)]}, 



where 



d/(^T)-[/K)^-g(^T)^]/(^±) 

D{u;) 



(4) 

(5) 
(6) 



resonance LUo{h) is simply ujQ^h) ^ loq ± h, where ujq c^l 
A(c7r/2)/ ln(8/7rc) is the resonance energy for h = O.Q 
Hence, for quasi- unitary scattering (wq <C A), already 
quite small values of h compared to A are sufficient to 
deeply modify the impurity-induced resonance. 

This is clearly seen in Fig.l where the w-dependence 
of No{lu) + (5Mmp(r, to) is plotted for h/A = 0, 0.02, 0.04. 
The value of the scattering parameter, c = 0.08, has been 
chosen in order to reproduce a zero-field maximum signal 
on the impurity site, r = (Fig. la), for w/^ = —0.03, 
i. e. the resonant energy reported in Ref.EI. For nonzero 
values of h, the resonance becomes Zeeman splitted in 
good agreement with ujQ^h) ~ i ^- This is also true 
for the LDOS signals away from r = 0, as it is shown in 
Fig.lb where the LDOS is plotted for r = (4/fcF,0) {i. 
e., along the direction of the gap maxima). 

The spatial dependence of the LDOS as a function of h 
for uj/A ~ —0.03 is shown in Fig.^ The pattern shown 
in Fig.||a {h = 0) closely resembles the spatial depen- 
dence obtained by Haas and Makij3 and it is character- 
istic of an electron-like bound state. However, -already 
for h/A = 0.02 (Fig.|b), which for A = 44 meVQ corre- 
sponds to a magnetic field of about 15 T, the resonance 
acquires a predominant hole-character, signalled by the 



Diu;) ^ {cl - [/(c^_) +5(c._)][/(c.+) + .9(c^+)]} 

x{cl - [/(c._) -g(c._)][/(c.+ ) -ffK)]} (7) 



where Cg, 



5(^±) 



l/iirNpSgVin-.p) = c/5g and 
d0 sin(0)^ 



27r [A((/.)2 



5c(r,w±) 
/c(r,'^±) 



d(f> cos(0) 



27r [A(0)2 



A(0) 



Aid,) 



(8) 
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The expressions for gs(r,Lu±) and fs(v,uj±) are obtained 
by replacing cos((/)) with sin(0) in the right-hand side of 
Eq.®. 

There are two important general features characteristic 
of the spin-orbit LDOS. First, as it can be inferred from 
Eqs.(yj9|), the spin-orbit LDOS vanishes at the impurity 
site: SN^{0,uj) = 0. Moreover, at /i = 0, SNsoir,0) = 
for every r. This is due to the k x k' factor appearing in 
Eq.(|l]) which makes the spin-orbit coupling particle- hole 
symmetry conserving. The second important feature is 
that the spin-orbit interaction can induce additional res- 
onances driven by the zeroes of Eq. . Without entering 
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FIG. 3. Effect of the spin-orbit contribution to the LDOS 
at r = (2/fcF,2/fcF) for c = 0.08 and different values of the 
spin-orbit parameter &g for /i = (a) and /i/A = 0.1 (b). 
Solid lines: &g = 0.01; dotted lines: ig = 0.05; dashed lines: 
&g = 0.1; long dashed lines: &g = 0.15; dot-dashed lines: 
&g — 0.2. The arrows indicate the positions of the coherence 
peaks for 5g = 0.01. 

too much into details, the main feature of the spin-orbit 
poles is thatj at ft- = 0, D{u!) vanishes at w = when 
Cso = l/7r.ll3 Away from this limit the poles acquire a 
finite imaginary part and move rapidly towards high en- 
ergies. Note however that since 5Nso{r,0) = 0, the reso- 
nance becomes sharper as Cgo —>■ 1/tt without reducing to 
a (5-function at = for Cgo = I/tt. Finally, due to the 
spin-mixing processes of the spin-orbit interaction, the 
effect of h is not merely a Zeeman split of the zero-field 
poles. 

These features are demonstrated in Fig.^ where the to- 
tal LDOS including the spin-orbit contribution is plotted 
as a function of u for r — {2/kp,2/kp), c = 0.08 and 
different values of dg. For h — (Fig.||a) and 6g — 0.01 
the presence of the coherence peaks at lo — ±A indi- 
cate that the LDOS is very close to that of a d-wave 
superconductor without impurities. However, as Sg is 
enhanced, the coherence peaks are depleted and a sym- 
metric broad resonance develops and moves towards low 
energies with contemporary reduction of its peak-width. 
The symmetry with respect to = merely refiects the 
particle-hole symmetry conservation of the spin-orbit in- 
teraction. At 6g — 0.2, the coherence peaks at w = ±A 
are completely suppressed and a sharp resonance is built 
at uj = ibwso with iOso ^ A. The origin of such low- lying 
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FIG. 4. Spin-orbit resonant energies lOso for c — 0.08 as a 
function of the external field. 

resonances stem from the poles of Eq.(0). Note in fact 
that for Sg — 0.2 the value of the spin-orbit scattering 
parameter, Cso = c/Sg = 0.4, nearly fulfils the condition 
Cso = I/tt for which, as discussed above, the spin-orbit 
T-matrix has a pole at LOgo = with vanishing imagi- 
nary part. The effect of the magnetic field is shown in 
Fig!|b, where the LDOS is plotted for h/A = 0.1 and 
for the same set of 6g values of Figj^a. As expected, for 
6g — 0.01 the coherence peaks at the gap edge are split- 
ted by ±/i. However, for higher values of Sg^ the low-lying 
spin-orbit resonances do not show Zeeman splitting be- 
cause of the presence of important spin-flip processes. 

Although the low-lying spin-orbit resonances are not 
Zeeman splitted, they nevertheless show some depen- 
dence on h. This is shown in Fig.^ where ujso is plotted as 
a function of h for Sg = 0.15 and 6g — 0.2. Note however 
that the /i-dependence of coso is rather weak, at least for 
low values of h. 

The spatial dependence of the total LDOS with Sg = 
0.2 is plotted in Fig.|| for the same parameters of Fig.^ 
(c = 0.08 and w/A = -0.03). For h ^ (Fig.|a) the 
poles of the spin-orbit T-matrix are at energies higher 
than Lu — —0.03 and the spatial dependence of the LDOS 
resembles closely that of Fig.^ where Sg = 0. Note 
however that the weight of the central peak is somewhat 
extended along the diagonals leading to an X-shaped ge- 
ometry. For h/A = 0.02 (Fig.||b) the spatial dependence 
is radically different from that shown in Fig.^. Now, a 
signal arises along the diagonals in the vicinity of r = 
with a contemporary shift of the peaks in the (±1, 0) and 
(0, ±1) directions at higher distances from the impurity 
site. This particular geometry of the LDOS is charac- 
teristic of the spin-orbit coupling to the impurity and, 
as inferred from Fig.^, it can be obtained also at h = 
when UJ is close to the spin-orbit T-matrix poles. 

In summary, two possible scenarios can be drawn about 
the Zeeman field effects on the LDOS of a (i^2_j^2-wave 
superconductor around a quasi-unitary impurity atom. 
First, if the spin-orbit coupling is absent or weak (cso ^ 
I/tt), the imposed magnetic field splits the quasiparticle 
resonance peaks by ±/i and, at fixed energy w, the spatial 
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FIG. 5. Spatial dependence of the LDOS with spin-orbit 
interaction for c = 0.08, Sg = 0.2 and uj/A — —0.03. (a): 
h = 0. (b): h/A = 0.02. 

dependence can be modified from electron-like to hole- 
like. Second, if the spin-orbit scattering is sufRciently 
strong (cso — I/tt), the LDOS acquires a novel off-site 
and particle-hole symmetric resonance at low energies 
which do not show Zeeman splitting at /i 7^ 0. Which of 
these two possibilities is actually realized would bring im- 
portant informations on the nature of the resonant states 
in high- Tc superconductors. 

A last remark regards the possibility of having very 
high (small) values of 6g (c) is such a way that Cso ^ 1- 
In fact, when Cgo — > 0, equations (ph and reduce 
to A^iij) = .g(co±)/[/(L^±)2 - 5(l^±P] and B±(^) = 
~f{^±)/[f{^±)^ ~ g{^±)'^\-, respectively. In this case, 
bN'^JyijS) does no longer contains spin-mixed terms and 
the two spin channels are perfectly decoupled. As ex- 
plained in Ref. |l^, this situation is due to the fact that, 
as long as the charge carriers are confined to move in 
the X — y plane, the spin-orbit impurity operator com- 
mutes with <Tz- For very strong spin-orbit interaction, it 
is found that also in the presence of an external mag- 
netic field perpendicular to the z direction, the (singlet) 
Cooper pairs are formed by electrons with opposite spins 
in the z-direction and the spin-rigidity of the supercon- 
ducting condensate is efficient against spin-flip transi- 



tions induced by the magnetic field. 

Note added - As shown in a recent publication [H. 
Tsuchiura, Y. Tanaka, M. Ogata, and S. Kashiwaya, 
Phys. Rev. Lett. 84, 3165 (2000)] a Zeeman spHtting 
of the quasiparticle resonances can be induced also by 
classical magnetic impurities. 
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